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ABSTRACT 

We study the stellar populations of a sample of 14 elliptical galaxies in the Virgo cluster 
Using spectra with high signal-to-noise ratio (S/N }t 100 A we propose an alternative ap- 
proach to the standard side-band method to measure equivalent widths (EWs). Our Boosted 
Median Continuum is shown to map the EWs more robustly than the side-band method, min- 
imising the effect from neighbouring absorption lines and reducing the uncertainty at a given 
signal to noise ratio. Our newly defined line strengths are more sucessful at disentangling the 
age-metallicity degeneracy. We concentrate on Balmer lines (H/3,H7,H(5), the G band and the 
combination [MgFe] as the main age and metallicity indicators. We go beyond the standard 
comparison of the observations with simple stellar populations (SSP) and consider four dif- 
ferent models to describe the star formation histories, either with a continuous star formation 
rate or with a mixture of two different SSPs. These models improve the estimates of the more 
physically meaningful mass-weighted ages. Composite models are found to give more con- 
sistent fits among individual line strengths and agree with an independent estimate using the 
spectral energy distribution. A combination of age and metallicity-sensitive spectral features 
allows us to constrain the average age and metallicity. For a Virgo sample of elliptical galax- 
ies our age and metallicity estimates correlate well with stellar mass or velocity dispersion, 
with a significant threshold around 5 x IO^^A/q above which galaxies are uniformly old and 
metal rich. This threshold is reminiscent of the one found by Kauffmann et al. in the gen- 
eral population of SDSS galaxies at a stellar mass 3 x IQ^'^Mq. In a more speculative way, 
our models suggest that it is formation epoch and not formation timescale what drives the 
Mass-Age relationship of elliptical galaxies. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: stellar content. 



1 INTRODUCTION 

Unveiling the star formation histories of elliptical galaxies is key 
to our understanding of galaxy formation. Being able to resolve 
their seemingly homogenous distribution is hampered by the fact 
that their light is dominated by old, i.e. low-mass stars, which 
do not evolve significantly even over cosmological times. Further- 
more, the presence of smal l amounts of young s t ars as recently 
discovered in NUV studies jperreras & SMbOOd : lYi et alJIlOOsI : 
lKaviraietal.]|2007h reveals a complex history of star formation 
that requires proper estimates of mass-weighted ages, in contrast 
with the luminosity-weighted ages that simple stellar populations 
(SSPs) can only achieve. The majority of papers dealing with age 
estimates of the stellar populations of elliptical galaxies rely on 
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such SSPs (see e.g. iKuntschner & Daviesll998l : [Trager et al.l200d : 
iThomas et al.ll2005h . and it is onlv recentlv that special empha- 
sis has been made on the need to go beyond simple popul ations 
jperreras & Yil2004l : ISerra & Trage3l200'7l : lldiart et alfcOOTl) 

Dating the (old) stellar populations of elliptical galaxies has 
been fraught with difficulties, the most prominent being the age- 
metallicity degeneracy, whereby the photo-spectroscopic proper- 
ties of a galaxy of a given age and metallicity can be replicated 
by a younger or older galaxy at a suitably higher or lower metallic- 
ity, respectively. To some extent, this problem has been overcome 
by the measurement of pairs of absorption line indices (* Worthevl 
1994), one whose change in equivalent width (EW) is dominated 
by the average metallicity of the population and the other domi- 
nated by the average age of the population. Typical metal-sensitive 
line strengths are the Mg feature at 5170 A , the iron line s around 
5300A, or a combination such as [MgFe] ( lGonzaleJl993h . Balmer 
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lines are more age-sensitive and are often combined with metal- 
sensitive lines to break the degeneracy. However, measurements of 
EWs of B aimer lines can be affected by the age-metallicity degen- 
eracy because of the presence of nearby absorption lines. Such is 
the case of H7, with the prominent G band at 4300A or the CN 
bands in the vicinity of the H(5 line. This paper is partly motivated 
by the need to define a method to estimate EWs that minimise the 
sensitivity of metallicity on Balmer lines by a proper estimate of 
the continuum. 

Even at relatively younger ages (a few Gyr) where such un- 
certainties are reduced, considerable degeneracies still remain. The 
confirma tion of recent star formation (RSF ) occuring in early type 
galaxies jYi et alj2005l : lKavirai et al.l2007t) has raised the problem 
that the existence of a young population can considerably affec t 
the para meters derived through SSP analvsis. iTrager et alj jlOOOh 
and later ISerra & TrageJ ( l2007f) showed that even relatively small 
mass fractions (~1%) of young stars can distort age and metallicity 
estimates and in moderate cases (~10%) completely overshadow 
the older population. In addition, the age and mass fraction of any 
younger sub-population will also be degenerate, with larger mass 
fraction of relatively older sub populations having the same effect 
as s maller fractions o f youn ger ones. 

ISchiavon et al.l ( |2004[) noticed that using different Balmer 
lines (H/3, H7 and H5) to estimate the age gives slightly different 
results, which was suggested to show tha t the galaxy ha d under gone 
recent star formation. Contrary to this, iThomas et al.l ( |2004|) find 
that H7 and H(5 equivalent widths are more affected by a non-solar 
a/Fe ratio on higher order Balmer lines. This effect is due to the 
increase of metal lines at bluer wavelengths, thereby distorting the 
co ntinuum as measured b y a side-band method. This was expanded 
byi erra & Tragej ( [20071) who remodelled synthetic 2-burst mod- 
els using H/9 and H7A and achieved a result consistent with both 
papers, concluding that a mistmatch between the three Balmer line 
estimates could possibly reveal underlying younger populations. 

Moving forward with such analysis - beyond simple popu- 
lations and luminosity-weighted parameters - requires improve- 
ments on both the H7 and YIS measurements. Balmer line equiv- 
alent widths suffer from the effects of the metallicity due to the 
presence of such lines in the spectral region used to determine 
the continuum ( Worthey & Ottaviani 1997; Thomas et al. 2 0041 : 
IProchaska et al.l l2007h . Althou gh considerable work has already 
been done in this area (see e.g. iRose et al.ll 19941: Jones & WorthevI 



1 19951 : IVazdekis & Arimotolll999l: lYamada et al.ll2006ir we test a 
completely different approach. 

In this paper we present a comprehensive analysis of the stel- 
lar populations of 14 Virgo cluster ellipticals using several models 
to describe the star formation history, exploring the discrepancies 
found between simple and composite stellar populations. In an at- 
tempt to combat the problems discussed above, we introduce a new 
method for the measurement of equivalent widths using a high per- 
centile running "median" to describe the continuum. The proper- 
ties of this new method are exploited by using various age- and 
metallicity-sensitive spectral features. 



2 THE SAMPLE 

We use a sample of 14 elliptical galaxies in the Virgo cluster, for 
which moderate resolution spect roscopy is available at high signal- 
to-noise ratio (S/N > 100 A~ MYamada et al.ll2006l . 120081 . Eight 
galaxies were observed with FOCAS at the 8m Subaru telescope: 
the other six were observed with ISIS at the 4.2m William Her- 



schel Telescope (WHT). Observations from Subaru span the spec- 
tral range A ~ 3800 — 5800A , whereas the spectra taken at the 
WHT span a narrower window, namely A ~ 4000 — 5500A . 
The resolution (FWHM) of both data sets is similar : 2A (Subaru) 
and 2 .4A(WHT). We refer the interested reader to lYamada et al.l 
( I2OO6I) for details about the data reduction process. We compare 
th ose spectra with compos ites of the R~ 2000 synthetic models 
of B ruzual & Ch ariot ( 2003), updated to the 2007 version ( Bruzuaj 
.20071) . We resampled the observed spectra from the orig inal 0.3 A 
to I A per pixel, performing an average of the spectra over a 1.5 A 
window, in order to have a sampling more consistent with the actual 
resolution. 

In this paper we use two alternative sets of information, either 
the full spectral energy distribution or targeted absorption lines. For 
the former, we consider a spectral window around the 4000A break, 
which is a strong age indicator (albeit with a significant degeneracy 
in metallicity, especially for evolved stellar populations). In order 
to minimise the effect of an error in the flux calibration, we do not 
choose the full spectral range of the spectra, restricting the analysis 
to 3800-4500A for the Subaru spectra, and 4000-4500A for the 
WHT spectra. 

The second method focuses on a reduced number of 
spectral lines. Following the traditional approach (see e.g . 



Kuntschner & Davies ll99i:lTrager et alj2000l : lThomas et alj200:. 
Sanchez-Blazquez et alj |2006h . we use a set of age-sensitive and 
metallicity-sensitive lines. In the next section we describe in detail 
the indices targeted by our analysis and describe a new algorithm 
that improves on the "standard" method to determine the contin- 
uum in galaxy spectra. 



3 MEASURING EQUIVALENT WIDTHS 

We focus on a reduced set of absorption lines originally defined in 
the Lick/IDS system ( Worthev et al. 1994) and extensions thereof 
dWorthev & Ottavianil [l997h . As age-sensitive lines we use the 
Balmer lines H/3, H7, H5, the G-band (G4300) and t he 4000A 
break (D4000). We use the standard definition of [MgFe] l lGonzale3 
Il993h . as a metal-sensitive tracer, which is a reliable proxy of over- 
all m etallicity, with a ver y mild dependence on [a/Fe] abundance 
ratio jThomas et al]|2003h . 

The standard method to determine the equivalent widths of 
galaxy spectra relies on the definition of a blue and a red side-band 
to determine the continuum. A linear fit to the average flux in the 
blue and red side-bands is used to track the continuum in the line 
(see e.g. iTrager et aljlioool) . This method, although easy to imple- 
ment, has an important drawback as neighbouring lines can make a 
significant contribution to the flux in the blue and red passbands, in- 
troducing unwa nted age/metallic ity effects. For instance, the H7 in- 
dex ( Worth ev & Ottavianill 19971) is defined with the blue side-band 
located close to the prominent G-band, around 4300A. This defini- 
tion causes non-physical negative values of the H7 line in absorp- 
tion, as the depression caused by the G-band makes the flux in the 
H-y line (wrongly) appear in emission. This has not prevented the 
community from using this line as a sensitive age tracer, as long as 
models and data are treated in the same way. IVazdekis & Arimotol 
( I1999I) defined new measurements of this line in order to reduce the 
metallicity degeneracy mainly introduced by the choice of the side 
bands. They avoid this by selecting specific regions less affected by 
the metal absorption lines. 

H(5 is another Balmer line which has been recently considered 
to be affected by neighbouring metal lines - most notably the CN 
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Figure 1. Dependence of the age-metallicity gradient - (A Age/Age)(/A Z/Z) - on the choice of BMC parameters, namely the level at which the 'boosted 
median' is taken (left) or the size of the kernel (right). On the left panel a number of kernel sizes is shown as labelled: AA=50A (dashed); lOOA (solid) and 
200A (dotted). On the right panel a number of confidence levels are considered: 50% (dotted); 70% (long dashed); 90% (solid) and 95% (short dashed). The grey 
horizontal line in both panels is the estimate from the standard side-band method (SB), and the grey vertical shaded area marks our choice of 'boosted median' 
parameters. 



molecular bands - wh ich reduce the age sensitivity of the index 
jProchaska et alj|2007h . 

In this paper we present an alternative method to determine 
the equivalent widths. Our method does not rely on the definition 
of blue and red side-bands and minimizes the contamination from 
neighbouring lines. This method is simple to apply and we propose 
it for future studies of stellar populations in galaxie|^ 

3.1 The Boosted Median Continuum (BMC) 

Our measure of equivalent width follows the standard procedure 
comparing observed flux in the line and the corresponding "inter- 
polated" continuum in the same wavelength range. For an equiva- 
lent width measured in A: 









L $c(A)J 



where Ai and A2 define the wavelength range of the line, "1>(A) 
is the observed flux, and <J>c(A) is the flux from the continuum. 
Rather than defining the continuum as a linear fit between a blue 
and a red side-band, we propose the "boosted median" of the flux, 
defined at each wavelength as the 90th percentile of the flux values 
within a lOOA window. 

This method is defined by two parameters, namely the choice 
of percentile (90% in our case) and the size of the kernel 
(AA =100A). The kernel size needs to be large enough to avoid the 
small scale variations of the spectra, but small enough to avoid dis- 
tortions from large scale structure of the spectra such as the breaks 



at 4OOOA and 4300A, or flux calibration errors 0. The choice of 
percentile also suffers a similar balancing act, since it should be 
high enough to select the true continuum but at a value that would 
avoid it becoming dominated by noise. In order to determine the 
optimal choice, a range of values for these two parameters was 
studied o n a number of simp l e stell ar populations taken from the 
models o f iBruzual & CharlotI ( I2OO3I) including the effect of veloc- 
ity dispersion and noise. Out of the simulations, we adopted the 
90th percentile of the flux within a lOOA window. However, this 
choice is not critical, given the robustness of the method in which 
the continuum is selected (i.e. a median). The advantage of averag- 
ing over a large enough wavelength range is that the effect of strong 
metallic lines in the vicinity of the index is limited. Furthermore, 
this pseudo-continuum is found to be less susceptible to noise (see 
below). 

Figure [T] motivates our choice of parameters. We show the 
age-metallicity sensitivity - (AAge/Age)/(AZ/Z) ~ for simple stel- 
lar populations measured at 12 Gyr and sola r metal licity. We do 
not follow the same definition as in IWorthevI ( fT994h . Instead, we 
change the age from the reference value by 2 Gyr and find the 
change in metallicity required from the fiducial SSP that gives the 
same variation in the EW. The grey horizontal bar is the value de- 
termined from the standard side-band method. Smaller values of 
the gradient imply a better disentanglement of the age-metallicity 
degeneracy. In the left (right) panels the horizontal axis explores a 
range of confidence levels (kernel sizes). The lines correspond to 
various choices of kernel size (left) and confidence level (right) as 
labelled. H/3 (top) behaves quite robustly with respect to the choice 



^ A C-programme that computes EWs using our proposed BMC method 
from an ASCII version of an SED can be obtained from us (fer- 
reras @star.ucl. ac.uk). 



^ Applying a flux calibration distortion of 5% - consistent with that found 
in current surveys (SDSS) - causes negligible effects on the EW measured 
with the BMC method (below 1%). 
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Table 1. Equivalent Widths of Virgo Elliptical galaxies using a 90% Boosted Median Continuum (see text for details). All values given in A and measured at 
the observed cr, with the Ict uncertanities in brackets below each measurement. 



Galaxy H/32o H720 H<52o Mgb2o Fe52702o Fe53352o G430020 D4000 [MgFe]2o 



NGC 4239 


82 


2.993 


2.097 


2.356 


2.972 


2.864 


1.673 


6.130 


1.441 


2.596 






(0.051) 


(0.064) 


CO 0561 


(0.044) 


(0.046) 


in 0601 


(0.046) 


(0.006) 


(0.032) 


NGC 4339 


142 


2.666 


1.488 


1.442 


3.807 


3.163 


1.772 


6.642 


1.560 


3.065 






(0.047) 


(0.054) 


(0.082) 


(0.045) 


(0.037) 


(0.050) 


(0.042) 


ro 0051 


(0.024) 


NGC 4365 


245 


2.191 


0.732 


1.195 


3.819 


2.971 


1.615 


6.289 




2.959 






(0.026) 


(0.042) 


(0.070) 


(0.023) 


(0.028) 


(0.031) 


(0.021) 




(0.015) 


NGC 4387 


105 


2.431 


1.698 


1.974 


3.614 


3.237 


1.703 


7.317 




2.988 






(0.044) 


(0.066) 


(0.068) 


(0.039) 


(0.057) 


(0.061) 


(0.050) 




(0.024) 


NGC 4458 


104 


2.464 


1.622 


1.676 


3.520 


2.719 


1.670 


6.570 


1.504 


2.779 






(0.039) 


(0.053) 


(0.055) 


(0.037) 


(0.042) 


(0.047) 


(0.041) 


(0.006) 


(0.025) 


NGC 4464 


135 


2.302 


1.547 


1.523 


3.648 


2.913 


1.521 


6.813 




2.844 






(0 0381 


fO 0531 


(0.070) 


(0.035) 


(0.034) 


(0.036) 


(0.037) 




(0.023) 


NGC 4467 


75 


2.605 


1.797 


1.875 


3.839 


3.209 


1.747 


6.760 


1.458 


3.085 






(0.050) 


(0.063) 


(0.065) 


(0.044) 


(0.052) 


(0.056) 


(0.049) 


(0.005) 


(0.026) 


NGC 4472 


306 


2.225 


0.793 


1.149 


3.650 


2.719 


1.535 


5.742 


1.525 


2.787 






(0.031) 


(0.069) 


(0.081) 


(0.039) 


(0.039) 


(0.042) 


(0.027) 


(0.005) 


(0.021) 


NGC 4473 


180 


2.410 


1.233 


1.592 


3.990 


3.212 


1.808 


6.918 




3.164 






(0.028) 


(0.044) 


(0.060) 


(0.026) 


(0.023) 


(0.041) 


(0.020) 


(— ) 


(0.021) 


NGC 4478 


132 


2.627 


1.653 


1.855 


3.719 


3.378 


1.838 


6.764 




3.114 






(0.042) 


(0.051) 


(0.069) 


(0.036) 


(0.039) 


(0.040) 


(0.039) 


(— ) 


(0.023) 


NGC 4489 


73 


3.204 


2.084 


2.328 


3.217 


3.248 


1.903 


6.455 


1.520 


2.878 






(0.047) 


(0.054) 


(0.075) 


(0.046) 


(0.048) 


(0.045) 


(0.050) 


(0.005) 


(0.028) 


NGC 4551 


105 


2.663 


1.717 


1.736 


3.907 


3.416 


2.027 


6.828 


1.531 


3.261 






(0.046) 


(0.058) 


(0.056) 


(0.034) 


(0.038) 


(0.044) 


(0.036) 


(0.005) 


(0.024) 


NGC 4621 


230 


2.284 


0.964 


1.312 


4.072 


3.014 


1.678 


6.250 




3.091 






(0.016) 


(0.024) 


(0.040) 


(0.012) 


(0.017) 


(0.018) 


(0.011) 


(— ) 


(0.009) 


NGC 4697 


168 


2.397 


1.289 


1.661 


3.806 


3.198 


1.832 


6.603 


1.524 


3.094 






(0.029) 


(0.037) 


(0.071) 


(0.027) 


(0.026) 


(0.033) 


(0.021) 


(0.004) 


(0.017) 



Velocity dispersions given in km/s, from lYamada et aj] j2006l) . 
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Figure 2. Comparison between the EWs measured by the side-band method 
(SB; horizontal) and our proposed 'Boosted Median Continuum' (BMC; 
vertical) for H/3, U-yp, H5p, G4300, Fe5270, Fe5335, Mgb, D4000 and 
[MgFe]. The error bars are shown at the Icr level. Notice only 8 galaxies 
(observed with Subaru) have a measurement of the 4000A break. For the 
remaining six galaxies our data does not extend bluer than 4000A. 



of BMC parameters. H7 (middle) shows that the presence of large 
scale features such as the break found around the G band at 4300A 
can affect the estimate if a large kernel size is chosen (200A, dotted 
line, left). Finally, 115 (bottom) shows that too small a kernel size 
(AA =50A, dashed line, left) or too high a confidence level (95%, 
short dashed line, right) can affect the estimate. In this case the ef- 
fect is caused by nearby features that will contaminate the estimate 
of the pseudo-continuum. Our choice of lOOA kernel size and 90% 
level is thereby justified by the need to avoid both short and long- 
scale features in the SEDs. Furthermore, lower confidence levels 
should be avoided as they will give flux values closer to an aver- 
age that will not reflect the true continuum given the presence of 
numerous absorption lines. Too high values of the level will make 
the measurement more prone to higher uncertainties at low signal- 
to-noise ratios. Monte Carlo simulations of noise show that for our 
choice of BMC parameters the errors for the BMC EWs are always 
smaller than those from the standard side-band method. 

Hence, the method for generating the BMC is fairly simple. 
At every wavelength the 90th percentile from all flux measurements 
within 50A on either side is assigned as the continuum at that wave- 
length. With the continuum thus defined, for each line strength we 
only have to define the central wavelength and the spectral win- 
dow over which the line is measured. Given that our method maps 
quite robustly the underlying continuum, we decided to fix a 20A 
width for all line strengths considered in this paper - hereafter, our 
BMC-based equivalent widths are labelled with a 20 subindex. 

This is motivated by the fact that for some definitions of the 
line strengths, any 'contaminating' lines falling within the central 
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Table 2. Range of parameters explored in this paper (see text for a description of each model). 



Model/Param 


MliN 


iVlAA 


Comments 


SSP 






2 params 


Age 


3 


13 


Gyr 




1 c 

-I. J 


+U. J 


Metallicity 


EXP 






3 params 


log T (Gyr) 


-1 


-1-0.9 


Exp. Timescale 


Zp 


0.1 


5 


Formation epoch 


log(Z/Z0) 


-1.5 


+0.3 


Metallicity 


2BST 






4 params 


to 


3 


13 


Old (Gyr) 


tY 


0.1 


3 


Young (Gyr) 


fY 


0.0 


0.5 


Mass fraction 


\og(Z/ZQ) 


-1.5 


-hO.3 


Metallicity 


CXP 






3 params 


logri (Gyr) 


-2 


-hO.5 


SF Timescale 


log T2 (Gyr) 


-2 


-hO.5 


Enrichment Timescale 


Zp 


0.1 


5 


Formation epoch 



bandpass can have a stronger effect on the BMC method compared 
to the standard side-band method. A clear example is found in the 
Hf3 line, where the effect of Fel( 4871 A) which sits on the shoul- 
der of the B aimer line within the standard definition of the central 
passband (width 28.75A ) clearly distorts the measurement of the 
equivalent width. This is slightly 'compensated' in the side-band 
method through the presence in the blue side-band of another sig- 
nificant Fe line at 489lA. This is clearly not the case for the BMC 
method, in which the pseudo-continuum is effectively independent 
of the values of both Fe absortion lines. Hence, in the BMC method 
it is desirable to choose a central passband which only targets the 
line of interest. Nevertheless, the method is versatile enough to de- 
fine wider central bandpasses. 

Thus, to avoid this problem we define a central pass- 
band wavelength as narrow as possible within the usual spec- 
tral resolutions targeted in unresolved stellar populations. We 
follow the 'F' type passb ands as used for the Balmer indices 
dWorthev & Ottavianilll99"7h . defining a 20A window centered on 
the line of interest for all indices consideretfE In the case of the 
metal line strengths - which usually involve clusters of lines - 
the major feature is chosen as the center of the new 20A index. 
In the case of Mgb2o, this is centered between the Mgl doublet 
at 5167/5172A. Fe52702o is defined by the Fel line at 5270A and 
Fe53352o uses the Fel line at 5328A. The central passband of the 
G-band is 4300A. Such a simple approach is nevertheless very ver- 
satile in its definition of any line st rength. The new D4 000 break 
feature uses the definition given bv lBalogh etai] ( Il999l) . with the 
difference that it is the ratio of the continuum flux obtained by the 
BMC method within those wavelengths that is used here. 

Table [T] shows the (BMC-measured) equivalent widths of 
the Virgo elliptical galaxies targeted in this paper. The measure- 
m ents are obtained dir ectly from the observed spectra, presented 
m lYamada et al] ( |2006|) and have not been corrected with respect 
to velocity dispersion. Analogously to the standard method, one 
could either correct the observed EWs for the effect of velocity dis- 
persion, or use the observed EWs in the modelling. For the latter, 
one must then measure the model EWs on spectra with the same 

3 Our BMC-based EWS are labelled with a '20' subindex, e.g. H/320 



resolution and velocity dispersion as the targeted galaxy. We fol- 
low this approach in the paper. We note at this point that the ef- 
fect of velocity dispersion on the EWs of Balmer lines measured 
with the BMC method is increased slightly since the method uses 
spectral information to better define the pseudo-continuum, which 
is destroyed by an increasing cr. This effect of course affects any 
other m ethod which measur es line strengths from unresolved spec- 
tra (e.g. lVazdekis & Arimoto 1999) . The numbers in brackets cor- 
respond to the 1-cr uncertainty, obtained from a Monte Carlo sim- 
ulation that generates 500 realizations of each SED, adding noise 
corresponding to the SNR of the observations. Notice that 6 of the 
galaxies do not have a measured D4000 as they were observed over 
a spectral range that does not include the blue passband used in the 
definition of D4000. 



3.2 Comparison with the side-band method 

Figure |2]shows a comparison of the line strengths measured on the 
same spectra using the side-band (horizontal axes) and the BMC 
methods (vertical axes). Notice the departure from a simple linear 
relationship, mainly caused by the different way neighbouring lines 
affect the estimate of the pseudo-continuum. In figures [3] and |4] we 
illustrate in more detail the difference in the measurement of the 
equivalent width of H7 and H(5, looking more closely at the effect 
of nearby lines on the indices. In order to eliminate the dependence 
of the EWs on velocity dispersion we classify the sample into two 
subsets, according to velocity dispersion. Figure [3] considers only 
those galaxies with a ^ 150 km/s and figure |4] focuses on galax- 
ies with G ^ 200 km/s. In both cases the galaxies are smoothed 
to the maximum velocity dispersion of the subsample, in order to 
make a consistent comparison. We note that by comparing spec- 
tra smoothed to some maximum a one would obtain a false sense 
of agreement between the methods, since the effect of smoothing 
removes information from the spectra. 

The left panels of figures |3] and |4] compare the EWs of the 
galaxies in each subset (open black circles), but we focus on two 
galaxies in each case (solid squares). The galaxies are chosen be- 
cause they have a similar value of the EW using one method, and 
a significantly different value using the other method. The spec- 
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Table 3. Ages and metallicities of Virgo elliptical galaxies according to the four models used in this paper. Error bars quoted at the 90% confidence level. 
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tral regions targeted in those galaxies is shown in the rightmost 
panels. The straight short line intersecting the spectra is the side- 
band pseudo-continuum and runs between the centres of the flank- 
ing side band regions The continuous line running along the top 
of the spectra is the BMC continuum. These are shown in grey 
and black for each galaxy, as labelled. The spectra are normalised 
to allow us to overp l ot bot h SEDs. In the case of H7 we follow 
IVazdekis & Arimotq ( Il999f) and normalise at the peak of 4365A. 
For H(5 we normalise to an average flux value of 1 within the wave- 
lengths of the blue side-band, to highlight the effect of CN absorp- 
tion towards the red side of the line. 

In figure [3] we show those galaxies with a ^ 150 km/s. 
In the top panels, we focus on NGC 4458 and NGC 4467, since 
these galaxies give a similar side-band H(5 EW but very differ- 
ent BMC-based measurements. Looking in detail at the SEDs 
on the right, the red side-band of NGC4467 is more affected by 
increased absorptio n, mainly caused by the CN features which, 
as reported in .Yamada et al.l (2006), are stronger in NGC4467 
(CNi=0.110 mag), than in NGC4458 (CNi=0.074 mag). Notice 
that this does not strongly affect the BMC pseudo-continuum. In 
addition, the H(5 absorption observed in NGC 4467 is possibly un- 
de restimated rela t ive to NGC 4458. The iron abundance reported 
by I Yamada eTai] ( I2OO6I) is higher in NGC 4467 by -0.1 dex, an 
effect that would be hidden by the normalisation over the iron lines 
around ~4070A. In the bottom panels of figure |3] we investigate 
the H7 lines in NGC 4387 and NGC 4467 which have a similar 
EW using the BMC method, whereas they differ by ~0.5 A in 
the side-band method. The SEDs reveal that these two galaxies do 
indeed have similar Balmer absorption, but the side-band method 
is 'tricked' when comparing the flux with the continuum (straight 
line). The increased G4300 absorption in NGC 4387 pushes the 
pseudo-continuum down, which reduces the EW value incorrectly. 
This is consistent with the measured absorption at this resolution 
of NGC4467 (G430020 = 6.5 A) compared to that in NGC4387 
(G430020 = 7.2 A). 

Shown in figure |4] is the subsample with higher velocity dis- 
persion, comparing galaxies at cr ^ 200km/s. The top panel shows 
the considerable discrepancy between the side band and BMC 
methods on measuring YiS. Since the over-abundance of CN and the 
resulting absorption tends to increase with the mass of the galaxy 



(see e.g. I Sanchez- Blazguez et al.l2003l : lTbroba et alj2009h , the in- 
clusion of the more massive NGC 4473 serves to show how sig- 
nificant the effect can be. The ~0.6A difference in H(5 given by 
the side band method can be seen to be generated by the increase 
in absorption in the red side-band region ~4120A, of NGC4473. 
We no te again that the CN absorption, identified by I Yamada et alj 
( I2OO6I) is higher in NGC 4473 (CNi=0.138 mag), than in NGC 
4478 (CNi =0.077 mag). The bottom panels display a more subtle 
situation in which the effect of the surrounding metal lines is to re- 
verse the relative values of the H7 EWs. The increased absorption, 
again mainly coming from the G4300 feature, forces the pseudo- 
continuum to a lower value, resulting in a small EW. The effect, 
while small, serves to indicate how sensitive the side band method 
can be to the surrounding spectral features. 

Therefore, figures [3] and |4] show that our proposed method is 
more resilient to the effects of neighbouring absorption lines. To 
further illustrate this point, figure [5] compares the age-metallicity 
degeneracy of typical Balmer {left) and metal lines {right) mea- 
sured either in the standard way or using the BMC pseudo- 
continuum. The shaded areas indicate the difference in the EWs 
with respect to a reference stellar population: on the left it is the 
population with the same age at solar metallicity, and on the right 
the reference is a population with the same metallicity and 10 Gyr 
old, both references marked by the vertical lines. Black and grey 
shading correspond to the BMC and the side-band methods, re- 
spectively. The panels on the left explore the age-sensitive Balmer 
indices for a range of metallicities and the panels on the right show 
metal indices for a range of ages. Ideally, a perfect observable 
would result in a horizontal line (i.e. zero metallicity dependence of 
an age-sensitive line and vice- versa). The models span a wide range 
of ages and metallicities as shown in the caption. This figure shows 
that BMC-based measurements of EWs are less subject to the age- 
metallicity degeneracy than the side-band methods. This result is 
especially dramatic for H7 and H(5, for which the age-metallicity 
degeneracy drops from AEW/A log(Z/Z0) = —4.2 (side-band) 
to -1.9 (BMC) in H7 or from -3.9 (side-band) to -1.7 (BMC) 
in H5 (values measured at the fiducial 10 Gyr, solar metallicity 
SSP). Furthermore, the shaded regions of the metal-line indices on 
the rightmost panels are much wider for the BMC method (black), 
showing that at a fixed age, the BMC method spans a wider range 
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Figure 3. Comparison of the H7i? and KSp equivalent width measured 
with the side-band (SB) and the Boosted Median method (BMC). In order to 
eliminate the dependence of the EWs on velocity dispersion, we select those 
galaxies with cr <^ 150 km/s and smooth them to this maximum velocity 
dispersion. Top: We select NGC 4458 & NGC 4467, which show similai" 
EWs in the SB method but have differing H5 EWs when using BMC. In the 
rightmost panels we plot the spectra from both galaxies. Overplotted are 
the pseudo-continuum of the side band (slanted straight line) and the BMC 
(running along the top of the SED). The identical EWs according to the SB 
method are due to the decreased flux in the red passband of NGC 4467, 
most likely caused by increased CN absorption, which hides the stronger 
intrinsic absorption in H<5. The spectra are normalised to an average value of 
1 across the blue passband oiYlSp to highlight the effect of CN absorption. 
Bottom: We highlight H7 in NGC 4387 & NGC 4467, for which the SB 
method indicates a considerable difference in EW while the BMC identifies 
almost none. Looking at the spectra (right) of both galaxies we can see that 
the difference in EW from the SB method is caused by an increase in the 
depth of the G4300 feature, artificially lowering the pseudo-continuum. The 
spectra are normalised to 1 at 4365A. 



Figure 4. As in figure[3]we compare the H'yp and HSp equivalent widths 
measured with the side-band and the Boosted Median method. In order to 
eliminate the dependence of the EWs and spectra on velocity dispersion, 
we select - in contrast to figure [3]- those galaxies with a 5C 200 km/s. 
The SEDs are smoothed to this maximum velocity dispersion. Top: The 
galaxies NGC 4473 & NGC 4478 are selected since they show differing 
H(5 SB EWs but similar BMC measured EWs. The spectra (right) identifies 
the cause of the discrepancy. The spectra of NGC 4473 redwards of the H<5 
fine is severely affected by CN absorption, lowering the average flux in this 
side band and distorting the derived pseudo-continuum. The BMC recovers 
a significantly more robust value. The spectra are normalised to an average 
value of 1 across the blue passband of HSp to highlight the effect of CN 
absorption. Bottom: NGC 4458 & NGC 4478 show discrepant SB EWs. 
This difference can be seen on the right, where the spectra of both galaxies 
are displayed. The BMC correctly identifies the larger value of H7 in NGC 
4478, which is hidden in the SB method due to the increased absorption 
of G4300 feature located in the blue pass band of H7^. The spectra ai'e 
normalised to 1 at 4365A. 



of EW, thereby being more sensitive to metallicity (see caption for 
details). 

In order to quantify the effect of an uncertainty in the veloc- 
ity dispersion on the EWs, we compare in figure |6] the change in 
EW between a fiducial value (a =100km/s) and a range of veloc- 
ity dispersions (horizontal axis). The SB (BMC) values are shown 
in the top (bottom) panels for a number of absorption lines, as la- 
belled. One can see that the effect on BMC-measured EWs can be 
corrected much in the same as with the standard SB method. In ad- 
dition, for most cases the correction is smaller for BMC estimates. 

Furthermore, we have compared the effect of noise on our pro- 
posed BMC pseudo-continuum and found that for a wide range 
of signal-to-noise ratios (from 10 to 100 per A ) the uncertainty 
in the EW of all lines is always ~ 0.3 dex smaller than those 
obtained with the side-band method. However due to the lower 
dynamical range, as the BMC is forced to stay above the spec- 
tra, the S/N requirements are si milar to other methods. We fol- 
low lCervantes & Vazdekisi feOOSl) and define the S/N at which the 
method can distinguish ±2.5 Gyr at 10 Gyr. For Hy and HS we 
get a S/N of ~ 50, and for H/J the required S/N is ~ 60. 

Figure [7] shows the equivalent widths of several lines for a 



grid of SSPs corresponding to a velocity dispersion of 100 km/s 
(black) or 300 km/s (grey). The measured EWs (without any cor- 
rection for velocity dispersion) are shown as black dots, along with 
a characteristic error bar which mostly comes from the systematics 
(the observed uncertainties are much smaller). We emphasize that 
our fitting method generates grids corresponding to the measured 
velocity dispersion for each galaxy. One can see that BMC EWs 
(left panels) appear less degenerate than the grids using a side-band 
method (right panels). 



4 MODELLING THE SFH OF ELLIPTICAL GALAXIES 

The properties of the unresolved stellar populations of our sample 
are constrained by comparing the targeted equivalent widths with 
four sets of generic models that describe the star formation history 
in terms of a reduced number of parameters. It is our goal to assess 
the consistency of different sets of models in fitting independently 
the different spectral lines targeted as well as the full SED. The 
majority of studies in the literature (see e.g. Kuntschner & Daviej 



1 19981 : iTrager et al.ll2000l : ICaldwell et all2003l : iThomas et alj2005l) 

have compared measurements of EWs with simple stellar popu 
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Figure 5. Dependence of Balmer lines on metallicity {left) and "metal" lines 
on age (right). The standard side-band method (SB; grey) and our proposed 
Boosted Median Continuum (BMC; black) are shown for a range of metal- 
licities and ages as shown. On the left, the shaded areas correspond to an 
age range [6, 12] Gyr. On the right, the shaded areas span a range of metal- 
licities: —0.3 < \o^Z/Zq < +0.3. The vertical axis is the difference 
between the equivalent width of the line for a given metallicity (left) or age 
(right) and the value at a reference point given by the vertical line (i.e. solar 
metallicity on the left and 10 Gyr on the right). 



Figure 7. Comparison of grids of SSP models using the standard side-band 
method (SB, right) and our proposed Boosted Median Continuum (BMC, 
left). A systematic error bar is included in each panel (the observational 
error bars are much smaller). We show the grids for two different velocity 
dispersions lOOkm/s (black) and 300km/s (grey). The solid lines connect 
the SSPs with ages (from top to bottom) 2 to 14 Gyr in steps of 2 Gyr. 
The dashed lines connect the SSPs with metallicities (from left to right) 
[m/H]=- 1 to +0.2 in steps of 0.2 dex. 
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Figure 6. The uncertainties caused by an eiTor in the estimate of velocity 
dispersion are illustrated by a comparison of the difference between the EW 
measured at a fiducial value (cr =100km/s) and a range of velocity disper- 
sions (horizontal axis). Various line strengths are considered, as labelled, 
for the standard SB method (top) and our proposed BMC method (bottom). 



lations (i.e. a single age and metallicity). While those models are 
probably valid for the populations found in globular clusters, it 
is imperative to go beyond simple stellar populations in galaxies, 
whose star formati on histories generate co mplex distributions of 
age and metallicity. iFerreras & Yj j2004h and lPasquali et al.l j2003h 
showed that composite models of stellar populations could re- 
sult in significant differences on the average ages and metallici- 
ties of galaxies. More recently, Sanch ez-Blazguez et al- ( 2006) and 
ISerra & Trager ( 2007) have explored this issue through the compar- 
ison of two age indicators, both concluding that composite popula- 
tions are needed to consistently model the populations in early-type 
galaxies. 

In this paper, except for the first case (namely Simple Stel- 
lar Populations), the models generate a distribution of ages and/or 
metal licities that are used to com bine the po pulation synthesi s mod- 
els of lBruzual & Chariot! Jiooj) , assuming a'Chabrier 72003') Initial 
Mass Function. We use the updated 2007 version ( Bruzual 2007i). 
The resulting synthetic spectral energy distribution is smoothed to 
the same resolution and velocity dispersion of the gal axy. A cor- 
rection for Galactic reddening is applied, assuming the lpitzpatriclj 
( Il999h law (this is mostly done for the comparison of the full SEDs, 
as EWs are not affected). Finally, the EWs are computed and com- 
pared with the observations using a standard maximum likelihood 
method. The high S/N of the observed spectra imply that our error 
budget is mainly dominated by the systematics of the population 
synthesis models. It is not obvious how to incorporate problems 
such as sparsely populated stellar parameter space or systematics 
within the models, into the analysis. As a compromise we esti- 
mate the uncertainties associated with the models by generating 
500 Monte Carlo simulations of gaussian noise at the level of S/N 
~50, that of the stellar library at the core of the BC03/CB07 mod- 
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Figure 8. Best fit age and metallicity values for our sample, using SSPs (left) or EXP models (right). The error bars shown the 68% confidence levels and the 
shaded regions give the age and metallicity estimates of Yamada et al. 1 2006) using the Hja vs. [MgFe] diagram. The diamond gives the fit to the age and 
metallicity using only the spectral energy distribution (i.e. no line strengths are used for this data point). 
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Figure 9. Best fit age and metallicity values for our sampl e, using 2-Burst (left) or CXP models (right). The en'or bars shown the 68% confidence levels and 
the shaded regions give the age and metallicity estimates of lYamada et al.l i200d) using the H7ct vs. [MgFe] diagram. The diamond gives the fit to the age and 
metallicity using only the specfi'al energy distribution (i.e. no line strengths are used for this data point). 



els ("STELIB lLe Borgne et al.l2003h . These are added in quadrature 
with those generated for each galaxy index. 

The four sets of models considered in this paper are listed be- 
low, with the range of parameters shown in table |2| 

(i) Simple Stellar Populations (SSP): This has been the most 
popular method used in the analysis of the ages and metallicities 
of early-type galaxies. The advantage lies in it simplicity: SSPs are 
the building blocks of all population synthesis models, and they are 
carefully calibrated against realistic SSPs (i.e. globular clusters). 



An SSP is primarily defined by an age and a metallicity although 
it is also dependent on the abundance pattern and the IMF. The 
drawback of this method is that the parameters obtained are inher- 
ently luminosity weiglited, such that a small amount of young stars 
can have a significant effect on the age extracted with this method. 
Furthermore, SSPs are not expected to model galaxy populations, 
which have an extended range of ages and metallicities. By using 
SSPs to model early-type galaxies one makes the assumption that 
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Figure 10. Likelihood distribution of the marginalised average stellar ages 
for some of the Virgo cluster elliptical galaxies from our sample. The dis- 
tribution is shown for the four models as labelled. The grey solid line cor- 
responds to the fit of the SED. The lines correspond to the age distribution 
according to [MgFe] plus either H/3 (black.solid), H7 (black,dashed) or H(5 
(black,dotted), respectively. 



the stellar populations have a narrow age distribution compared to 
stellar evolution timescales. 

(ii) Exponential SFH (EXP): A more physical scenario should 
consider an extended period of star formation. The EXP models 
(also called r models in the literature) model the star formation 
rate as an exponentially decaying function of time with timescale 
r, started at an epoch given by a formation redshift zp- The metal- 
licity is assumed to be fixed throughout the SFH and is also left as 
a free parameter. 

(iii) Two Burst Formation (2BST); Recent rest-frame NUV ob- 
servations have revealed the p resence of residual star formation 



in early-type galaxie s (see e.g. Trager et"^l2000l : iFerreras & Silkl 
I2OOOI : lYi et al.ll2003 : iKavirai et alj|2007h . The presence of small 
amounts of young stars can significantly affect the derived SSP 
model parameters tSerra & Tr ager 2007). This model describes this 
mechanism with two simple stellar populations. Four parameters 
are left free: the age of the old (to) and the young components 
(ty), the mass fraction in young stars (fy), and the metallicity of 
the populations {Z, assumed to be the same for the old and the 
young components). 

(iv) Chemically Enriched Exponential (CXP): All the models 
considered above fix the metallicity throughout the SFH. More 
realistically, a model should incorporate in a consistent way the 
buildup of metallicity caused by previous generations of stars. 
These chemical enrichmen t models (see e.g. iBoissier & Prantzoi 
ll999l : lFerreras & Silkl2003l) include the stellar yields from interme- 
diate and massive stars and result in a distribution of metallicities 
as well as ages. The aim of this paper is to explore simple mod- 
els that can be easily implemented. Hence, instead of applying a 
detailed model of chemical enrichment, we define in a purely phe- 
nomenological way a model that mimics those. We keep the same 
SFH as in the EXP models (described by a star formation timescale 
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Figure 11. Compaiison between the ages presente d here and tho s e from 
previous _studies. Y06, TF02 and C03 indicate 1' Y amada et all l2006l) . 
iTerlevich & Forbes 200^, and (Caldwell et al. 2003), respectively. The 
figure illustrates the difference in the average ages when including two types 
of composite populations: a 2-burst model (2BST; top) or a smooth star for- 
mation history including chemical enrichment (CXP; bottom). The slanted 
solid line is a 1 : 1 correspondence. The vertical dashed line represents the 
age of the Universe for our cosmology, and is used as a prior in the analysis 
(see text for details). 



ri and a formation redshift zp). Furthermore, the metallicity is as- 
sumed to increase in step with the cumulative distribution of the 
star formation rate, namely: 



Z{t) = + {Z2 - Zi] 



1 - exp(-At/r2) 



(2) 



where Ai = t — t{zF) and T2 is the timescale corresponding to 
the buildup of metallicity. As a first-order approach it is valid to 
assume that metallicity increases with the stellar mass of the sys- 
tem (roughly, the cumulative distribution of the star formation rate). 
More accurately, this timescale will depend on the star formation 
efficiency or on the fraction of gas ejected in outflows (see e.g. 
IFerreras & Sii3l2003h . The upper and lower values of metallicity 
are fixed to Zi = Zq/1G and Z2 ~ ^Zq, although the average 
metallicity will be controlled by the timescale T2. 

We assume a standard ACDM cosmology (£7^=0.3; Ho= 70 
km/s/Mpc) to constrain the maximum ages of the stellar popula- 
tions to the age of the Universe (13.7 Gyr). One could allow the fit- 
ting algorithm to stray into older ages to explore wider volumes of 
parameter space, but this paper takes the cosmological cons traints 
on the age of the Universe at face value dSpergel et al.ll2007h . Fur- 
thermore, population synthesis models are poorly calibrated for 
populations older than Galactic gl obular clusters. We note that pre- 
vious work on the same spectra ( Yamada e t al. 2006) probed the 
population synthesis models of IVazdekis & Arimoto. (,1999) all the 
w ay to their old e st ava ilable ages (~20 Gyr), using the isochrones 
of iBertelli et al.l ( Il994l) . Although in only two cases, NGC4365 at 
20 Gyr and NGC4464 at 18.5 Gyr, did they find cosmologically 
contradictory ages, we still choose to avoid this by limiting the 
maximum age. 



4.1 Disentangling the stellar populations 

The comparison of the EWs is done by measuring the line strengths 
directly from the model SEDs after being resampled and smoothed 
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Table 4. Best fit parameters for the 2BST and CXP models 
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1 Stellar masses are computed from the best fit CXP (2BST) models. Estimated uncertainty A logMa ~ 0.3 dex (mostly from the assumption about the IMF). 



to the resolution and veloci ty dispersion o f the actual observations 
of each galaxy, taken from Yamad a et al ] ll200^. Since the spec- 
tra of the models have already been smoothed in terms of velocity 
dispersion to that of the real galaxy, no correction for this need 
be applied. All galaxies are fitted using the age-sensitive indices, 
H/320, H720, H(52o and G4300 20 along with the metallicity indica- 
tor [MgFe] 20 ( Gonzalez" ' 1 993') . This is used instead of the newer 
[MgFe]' as defined by Thomas et all since the original is 

almost independent of a/Fe and it is not clear whether the same 
fractions of the two Fe lines is still appropriate for the new BMC 
indices. We note that we have excluded the D4000 index since 
it was found to be poorly modelled by B C03, a problem already 
identifi ed by both Wild et al. (JOQZ) and ISanchez-Blazquez et al.l 
( l2006ch . The exclusion (or inclusion) of the D4000, does not sig- 
nificantly affect the final results obtained although it does signifi- 
cantly increase (decrease) the quality of the fit. The removal of the 
D4000 also limits the constraints we can put on second order pa- 
rameters (e.g. age and mass fraction of young stellar population in 
the 2-Burst model). 



As indicat e d bylSchiavon et al. U2004h and later confirmed by 
ISerra & Trageil j2007h . a discrepancy found between the parame- 
ters derived through each of the three Balmer lines could be ex- 
plained by the existen ce of an underlying younger sub-population. 
ISchiavon et al.l j2004h explain that this may well be due to the dom- 
inance of the younger population at bluer wavelengths. However 
the mere fact that different indices would give different parame- 
ters, should suggest that the model being used has not captured 
all aspects of the true star formation history. Hence, in the anal- 
ysis of EWs of our sample galaxies, we are not only looking to 
determine the age and metallicity but also estimating whether it 
is possible to rule out a single population to describe an early- 
type galaxy. This could be done through discrepancies between 
the three Balmer line s targeted here, a possibility highlighted by 
ISchiavon et al. 1/2007) on a single stacked spectra of low mass el- 
lipticals, who found that the addition of a small young population 
alleviated such discrepancies. 

With the same aim, we focus on possible discrepancies be- 
tween individual line strengths and explore a number of star for- 



mation histories to determine the age distribution of our sample of 
elliptical galaxies. 



4.2 Spectral Fitting 

As an alternative approach to targeted line strengths, we also con- 
sider the full spectral energy distribution to constrain the star for- 
mation history of our sample. In this paper we want to test the con- 
sistency between the constraints imposed by the line strengths and 
those from a full spectral fit. 

We use a maximum likelihood method to fit the spectra, which 
is analysed over the spectral range between 3900A and 4500A (or 
4000-4500A for the six galaxies observed with WHT). The SEDs 
are normalised over the same range. The data span a much wider 
spectral range (out to 5500-5800A). However, we have chosen a 
smaller range to avoid flux calibration errors, which could intro- 
duce important systematic changes in the predicted ages and metal- 
licities. The spectral range chosen straddles the D4000 break, since 
this is a very sensitive indicator of the ages of the stellar popula- 
tions. 



4.3 Results 

Figures [8] and [9] show the average age and metallicity for all four 
models explored in this paper. The error bars are shown at the 69% 
confidence level (i.e. like a \a level although our analysis does not 
assume a Gaussian distribution). For the SSPs the average age is 
trivially the age of the population. For the others, the average val- 
ues are weighed by the stellar mass. Each error bar corresponds to 
the individual fit of an age-sensitive line (H/3, H7, H5 or G4300) 
plus the [MgFe] as metallicity indicator (i.e. each point comes from 
the analysis involving 2 indices). The diamond and its error bar cor- 
respond to the fit to the SED - no additional information from the 
line strengths is added to this likelihood. In grey, the shaded areas 
correspond to the ages from Yamada et al. ( 2006) (using the YL^fcr - 
[MgFe] diagram on SSPs). One can see that SSPs (left panel of fig- 
ure [8} show the largest discrepancies between the ages and metal- 
licities estimated using independently the different age-sensitive in- 
dices. This would imply that a joint likelihood putting all indices 
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together would not be a consistent way of determining tlie ages of 
the populations. 

Nevertheless, as a comparison we show in table[3]the best esti- 
mates of the average age and metallicity, combining all five indices 
(but not the fit of the SED). Table |4] shows the best fit parameters 
for the 2BST and CXP models. The uncertaintes in both tables are 
given at the 90% confidence level. It can be seen from table[3]that in 
a majority of cases the use of a model with an age spread offers an 
improvement in the value for the five indices. Importantly the 
recovered values of metallicity remain consistent (within errors) 
across the single metallicity models considered, indicating its ro- 
bustness against the inclusion of additional populations. For some 
galaxies, such as NGC 4365, NGC 4387 and NGC 4473, an SSP 
is as good as the other models presented here and in fact the pa- 
rameters recovered by complex SFHs essentially reduce to an SSP. 
This is to be expected since the short formation timescale expected 
for elliptical galaxies will in some cases be well represented by 
an SSP. However galaxies which have younger SSP ages give very 
different (mass-weighted) ages when considering a 2-Burst model. 
One interesting result is that the CXP models, which incorporate a 
spread both in metallicity and age, find solutions that are as good 
as - and in the case of NGC4458, and NGC4467 and NGC 4697 
significantly better than - the other models. 

In addition, the discrepancies among individual measurements 
shown in figure [8] reflects the fact that a single age and metallicity 
scenario is a weaker, less consistent model of the populations of a 
galaxy. Table[3]shows that the actual change in the age when going 
beyond SSPs can be quite significant, especially for younger ages. 
We emphasize here that the different ages do not reflect a possible 
inclusion of a prior caused by the choice of parameters. All the 
composite models presented here (CSP, EXP and 2-Burst) explore a 
range of parameters that includes the best-fit ages and metallicities 
obtained by the SSPs. 

In some cases the values are larger than might be expected. 
In the case of NGC 4621, the largest contributor is the G-band. 
If we remove that indicator from the analysis, the reduced de- 
creases to ~ 1. It is clear that the high S/N of the spectra poses a 
challenge for current stellar population synthesis models. Certainly 
the errors associated with the models do not take into considera- 
tion internal problems with the models themselves or the possible 
limitations of the stellar libraries from which they are assembled. 
IProctor et al. et al] ( l2004h also consider that limitations in the mod- 
els as the cause of the high values. 

The mismatch discussed above can be seen in more detail in 
figure [Tol where the marginalised distribution of average age for 
all four models is shown for a few galaxies from the sample, as la- 
belled. The black solid, dashed, and dotted lines correspond to in- 
dividual fits to H/3, H7, and H(5, respectively, along with the G4300 
(grey, dotted) and the SED (grey, solid). SSPs (leftmost panels) fail 
to give a consistent distribution, whereas any of the other models 
give a more unified distribution. Notice that the galaxies shown in 
the top panels (NGC 4239 and NGC 4489) have young populations, 
whereas the bottom two galaxies (NGC 4464 and NGC4697) have 
overall older populations. The top two galaxies are better fit by a 
2-Burst model and the bottom two get better fits from an extended 
model like CXP. This would suggest that the youngest populations 
in early-type galaxies are best fit by the assumption of "sprinkles" 
of young stars, as s ugges ted by Trager et al. (2000) and modelled 
bv lFerreras & Silkl bOOd) and lKavirai et alj ( l2007h . 

This explains why the models incorporating an age spread 
have greater maximum likelihood values, where Cmax oc 
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Figure 12. Ages and metallicities versus stellar mass (left) and velocity dis- 
persion (right). The solid dots (grey diamonds) correspond to the estimates 
from the CXP (2-Burst) models. The dashed lines in the ri ghtmost panels 
are the age and metallicity scaling relation of lThomas et al.ll2005i) . 



SSPs are not substantially worse than the more complex models. 
This is not surprising given firstly, that the SFH of many elliptical 
galaxies is not different from a single episode of star formation, 
es pecially when it happen s at high redshift. Secondly, as identified 
bv lSerra&Tra "^( l2007h . for bursts of mass fraction, /y ^ 10% 
the disagreement in the recovered ages from different Balmer lines 
dissappears. This indicates that there is no differential effect across 
the spectrum from which to determine an age spread. Nevertheless, 
it remains the case that the more complex models are physically 
well motivated, resulting in better average age estimates and fit the 
data more consistently. Under such criteria we select the 2BST and 
CXP for further analysis. 

Figure [TT] shows a comparison of our age estimates using 2- 
Burst (top) and CXP models (bottom) with values taken from the 
literature, as given in the figure caption. It is important to no- 
tice that our models do not consider ages older than the age of 
the Universe using a concordance ACDM cosmology (i.e. 13.7 
Gyr). Our ages are also closer to a mass weighted average due 
to the underlying assumptions in the modelling in both cases. 
We find that for the youngest galaxies, our analysis gives rather 
older ages than previous estimates from the literature. Also no- 
tice that the older galaxies appear younger because of the strong 
constraint imposed by the chosen cosmological parameters. Fig- 
ure [12] shows the average age and metallicity as a function of stel- 
lar mass ( left) and velocity dispersion (right), both for the CXP 
(black dots) and for the 2-Burst models (grey diamonds). The stel- 
lar masses are obtained by combining the absolute magnitude of 
the galaxies in the V band with M/Ly corresponding to the best 
fit CXP models (There is no significant difference if the best fit 
2-Burst models are used instead). The dashed line is the fit to 
age and metallicity from iThomas et all ( |2005|) . Our age-mass re- 
latio nship is compatible with recent estimates in the literature (see 
e.g. iTrager et alj|2000l; JCaldweU et alj|2003l ; iThomas et al]|2005l: 



exp(— Xmin/2). However, as was seen in table [3] the fits of the ISanchez-BIazguez et alj 2006bl) . We should emphasize that our re 
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Figure 13. The predictions of the CXP models for star formation timescale 
( top) and formation redshift (bottom) are shown with respect to stellar mass 
(left) and velocity dispersion (right). This speculative plot suggests that the 
age trend found in figure flll is mainly caused by a range of formation epochs 
(bottom). The formation timescale (top) also coiTelates with mass or ve- 
locity dispersion, but it does not have values larger than ri ^ 1.5 Guyr. 
The solid lines in the right panels are simple linear fits to the data (equa- 
tions[3]&|4] see text for details). 



duced sample of Virgo elliptical g alaxies excludes mas sive, young 
SOs in high density environments jXhomas et ai]|2005l) , removing 
them as a possible cause of the age scatter. 

Furthermore, notice in figure [12] the transition in the age dis- 
tribution around 150 km/s in velocity dispersion or 5x IO^^A/q in 
stellar mass. Above this value the ages of the populations are over- 
all old and with small scatter among galaxies, whereas lower mass 
galaxies have a wide range of ages and metallicities (not necessar- 
ily younger throughout). The trend is robustly independent of the 
modelling, as shown both by the CXP fits (black dots) and the 2- 
Burst models (grey diamo nds). This t hreshold is reminiscent of the 
one found by Kauffmann et al.l ( 12003 ) in the general population of 
SDSS galaxies at a stellar mass 3 x 10^^'' Mq or the threshold in 
star formation efficiency of late-type galax ies at rotation veloci ties 
Vc ~140 km/s based on their photometry jperreras et al.ll2004h or 
on the presence of dust lanes jPalcanton et alj2004l) . 



5 CONCLUSIONS 

We have revisited the superb, high signal-to-noise spectra of 14 
elliptic al gal axies in the Virgo cluster presented by Yamada et al. 
I2OO6LI2OO8I) . Our main goal was two-fold: a) to give an optimal 
but versatile definition of the equivalent width of spectral features 
that would minimise the contribution from nearby lines, b) to ex- 
plore the possibility of discriminating between the standard treat- 
ment of galaxy spectra either as simple stellar populations or com- 
posite models with a distribution of ages and metallicities. 

The former was achieved by defining a "boosted median" 
pseudo-continuum. This method is very easy to implement on any 
spectral data and it improves on the traditional side-band meth- 



ods by reducing the final uncertainty of the EWs for the same 
spectra and reducing the age-metallicity degeneracy of age- and 
metal-sensitive lines by reducing the effect of unwanted, nearby 
spectral features. The method only requires two numbers to deter- 
mine the pseudo-continuum (i.e. confidence level of the boosted 
median and size of the kernel over which the analysis is per- 
formed). We propose 90% and lOOA for these two parameters when 
dealing with stellar populations of galaxies at moderate resolution 
(R~ 1000 - 2000). 

The second goal - going beyond SSPs - is approached by 
comparing SSPs and three more sets of models which assume some 
distribution of ages and metallicities. We find simple populations 
fail to give a consistent marginalised distribution of ages when fit- 
ting independently different age-sensitive lines. We propose either 
a 2-Burst model or a r star formatiom history with a phenomeno- 
logical prescription for chemical enrichment (CSP models). They 
give similar results, with a clear trend between average age and ei- 
ther stellar mass or velocity dispersion, as shown in figure [72l 

As suggested above when discussing figure [TO] if younger 
(older) galaxies are better fit by 2-Burst (CXP) models, one would 
expect that young stars in elliptical galaxies appear in a random 
way, and not as a time-coherent, smooth distribution. This re- 
sult supports minor merging - possibly involving metal-poor sub- 
systems - as the main mechanism to generate recent star formation 



in early-type galaxies jKavirai et al.ll2009l) . This process will be 
readily detected in low-mass galaxies, whereas a similar amount of 
young stars in a massive galaxy will be harder to detect. 

Doubtlessly average ages and metallicities are the observ- 
ables best constrained by these models. Parameters like formation 
timescale and formation epoch can only be considered "next-order" 
observables, which will be fraught with larger uncertainties. This 
is largely due to the massive degeneracies which are present in 
such parameters as mentioned in the introduction. Its is well known 
that the mass and age of a subpopulation are degenerate in that 
additi onal mass fraction can compensate for an older age or vice 
versa i lLeonardi& Rose! [1996). Another considerable degeneracy 
also exists between the formation redshift and the st ar formation 
times cale for both the EXP and CXP models (see e.g. iGobat et"ail 
|2008f). While our indices cover a large wavelength range we find 
that these degeneracies are still significant. FigurefTSjshows the de- 
generacy between parameters for a typical case using 2-Burst mod- 
els (left) and CXP models {right). We take the probability weighted 
values across the entire parameter space as a more robust mea- 
sure of these parameters. Although the derived parameters are still 
subject to considerable uncertainty, as expressed by the error bars 
shown in figure [T3] we note that our grids sample the parameter 
space quite finely meaning our recovered likelihood distribution 
should be a reasonable approximation to the complete distribution. 

We proceed in a more speculative way, taking the continuous 
CXP models at face value and putting some trust on the predic- 
tions of their formation epochs and timescales. Figure [T3] reveals 
an intriguing trend that suggests it is not just formation timescale 
but also formation epoch which drives the mass-age relationship in 
early-type galaxies. The solid lines in the top/bottom right panels 
gives a simple linear fit to the data, namely: 



zf ~ 1.8 + 3.9 log (Jioo 
r(Gyr) ~ 0.3 — 1.14 log ctioo 



(3) 
(4) 



with (Jioo given in units of 100 km/s. Figure [14] shows our pro- 
posed model for the star formation history of Virgo cluster ellipti- 
cal galaxies. We model the formation histories as skewed Gaussians 
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Figure 14. Our proposed model for the star formation history of Virgo cluster elliptical galaxies. The parameter estimates suggest a nearly 
constant star formation timescale with respect to velocity dispersion, whereas formation epoch is strongly con'elated. We find a stronger 
correlation with velocity dispersion (compared to stellar mass), so we label the formation scenarios according to cr. The bottom panel 
emphasizes the issue that the obsen'ed low-mass Es are best fit by later formation epochs, an issue which does not rule out the possibility 
of other low-mass Es in the past - or elsewhere - to have been formed earlier (grey Gaussians). 



whose spread maps the star formation timescale, and we include a 
correlation between formation epoch and velocity dispersion, as la- 
belled. Notice that the observed low-mass galaxies (black in the 
bottom panel) are best fit by late formation epochs, although this 
fact does not prevent smaller early-type galaxies in the past to have 
formed at earlier epochs (represented by the grey Gaussians). Our 
sample is too small and too local to explore this issue further. 

One could argue that the lack of a correlation between for- 
mation timescale (n; top) and mass would contradict the ob- 
served relationship between mass and abundance ratios (see e.g. 
I Tho mas et al. 2005). One possibility would involve minor merging 
of dwarf galaxies. Those mergers will have a more prominent effect 
on the spectra of low-mass galaxies. Dwarf galaxies have a very ex- 
tended period of star formation and eject a big fraction of their met- 
als, resulting in metal-poor gas with solar or sub-solar abundance 
ratios that will reduce the observed (luminosity-weighted) [a/Fe] 
in low-mass galaxies. Alternatively, notice that the formation epoch 
(zf; bottom) correlates quite strongly with velocity dispersion, im- 
plying that the structures that form the bulk of the stars in low- 
mass galaxies are more likely to be contaminated by the ejecta from 



type la supernova, thereby reducing the abundance ratios to repro- 
duce the observed correlation between [a/Fe] and mass. A more 
detailed analysis of the abundance ratios - although beyond the 
scope of this paper - is under way to explore this very interesting 
possibility. 
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